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In many applications a source of the blackbody radiation (BBR) can be highly anisotropic. This 
leads to the BBR shift that depends on tensor polarizability and on the projection of the total 
angular momentum of ions and atoms in a trap. We derived a formula for the anisotropic BBR shift 
and performed numerical calculations of this effect for Ca^ and Yb^ transitions of experimental 
interest. These ions were used for a design of high-precision atomic clocks, fundamental physics tests 
such as the search for the Lorentz invariance violation and space-time variation of the fundamental 
constants, and quantum information. Anisotropic BBR shift may be one of the major systematic 
effect in these experiments. 

PACS numbers: 32.10.Dk, 06.30.Ft, 31.15.Am 


I. INTRODUCTION 

The past five years brought remarkable improvements 
in both accuracy and stability of atomic clocks m- 
Development of ultra-precise clocks is important for a 
wide range of applications, including design of abso¬ 
lute gravimeters and gravity gradiometers for geophys¬ 
ical monitoring and research, gravity aided navigation, 
improved timekeeping and synchronization capabilities, 
tests of fundamental physics such as Einsteins theory of 
relativity, search for variation of fundamental constants 
through time, space, or coupling to gravitational fields, 
and exploration of strongly correlated quantum many- 
body systems [5]. 

Performing fundamental tests with the atomic clocks 
and other precision atomic, molecular and optical (AMO) 
technologies leads to ever increasing requirements for 
the understanding and control of the systematic errors. 
Moreover, a number of novel fundamental physics AMO 
experiments including searches for ultralight (sub-eV) ax- 
ions, axion-like pseudoscalar and scalar dark matter [7HS] 
and topological defect dark matter HO] have been car¬ 
ried out or proposed, requiring improved understanding 
of systematic effects in AMO systems. 

One of the major experimental and theoretical prob¬ 
lems in improving the atomic clock accuracy is a preci¬ 
sion determination of atomic clock frequency shift due 
to blackbody radiation (BBR). In recent experimental 
works m with ®’^Sr optical lattice clocks, the blackbody 
radiation was identified as the primary source of clock’s 
uncertainties. A number of measurements and thorough 
analysis of all systematic effects led to reduction of the 
Sr clock total uncertainty to the level of 2.1 x 10“^® in 
fractional frequency units. However, 65% of the clock 
uncertainty budget was still due to the BBR shift [T]. 

A number of other experiments with trapped ions and 
atoms are sensitive to the BBR effects, including recent 
tests of local Lorentz invariance (LLI) violation in the 
electron-photon sector with trapped Ca+ [TT] ions. Theo¬ 


ries aimed at unifying gravity with quantum physics sug¬ 
gest that Nature violates Lorentz symmetry at the Planck 
scale while suppressing its violation at experimentally 
achievable energy scales m- The minimal 0(1) sup¬ 
pression may lead to Lorentz-violating effects appearing 
beyond 10“^^ sensitivity level, determined by the ratio 
of electroweak and Planck scales. Thus, high-precision 
experiments with atomic systems mi [13] provide an im¬ 
portant route to search for Lorentz violation at low en¬ 
ergies. The LLI experiments with trapped ions may be 
particulary sensitive to anisotropic BBR shift since they 
are based on monitoring the energy difference between 
different Zeeman substates as described below. In these 
experiments, anisotropic BBR shift may become a lim¬ 
iting factor for the ultimate accuracy of the Lorentz vi¬ 
olation tests in the electron-photon sector mi and this 
work is strongly motivated by these fundamental stud¬ 
ies. BBR effects may also become a source of deco¬ 
herence in larger-scale quantum information experiments 
with trapped ions due to a change in the environmental 
temperature or temperature gradients during the com¬ 
putation. 

Calculations of blackbody radiation shifts are usually 
done assuming that the BBR radiation is isotropic. 

A detailed consideration of the isotropic BBR effect 
in conventional electric dipole approximation was car¬ 
ried out in mi. Multipolar theory of isotropic BBR shift 
of atomic energy levels (as well as its implications for 
optical lattice clocks) was developed in [T3j. However, in 
practice the source of BBR can be highly anisotropic and 
even may have a small angular size. As a result, a lot of 
experimental efforts is required to make the BBR field 
uniform and with a known temperature m m [H] ■ For 
this reason a proper calculation of the anisotropic BBR 
shift is necessary. Also, this effect can be of interest by 
itself as a physical phenomenon. We note that the prob¬ 
lem of anisotropic BBR effect, discussed in the present 
work, does not arise in experiments with alkaline-earth 
atoms aiming to create an atomic clock in the 
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transition because the states with total angular momenta 
J = 0 are involved. 


II. GENERAL FORMALISM 

The isotropic BBR shift of an energy level is propor¬ 
tional to scalar static polarizability of this level [HHS]. 
In the case of anisotropic blackbody radiation, an addi¬ 
tional contribution that depends on the projection of a 
thermal photon wave vector k to the z axis arises. We 
show below that it is determined by the tensor polariz¬ 
ability of the level. This contribution is particularly im¬ 
portant when we consider BBR frequency shift of a \JM) 
- \JM') transition, where M and M' are the projections 
of the total angular momentum J to the z axis, and we 
assume J and M to be good quantum numbers for the 
atomic states. The scalar polarizabilities of the Zeeman 
substates \ JM) and \JM') are practically identical dif¬ 
fering only due to a very small difference in the energy 
denominators. As a result, the isotropic BBR frequency 
shift is completely negligible in this case. In contrast, 
the anisotropic BBR shift of the energy level depends on 
and can be noticeably different for the different M 
substates. The same issue arises for the hyperfine states 
with different Mp. 

Such a systematic effect arose in a recent record- 
high precision experiment aimed at the search for local 
Lorentz invariance violation in the electron-photon sec¬ 
tors using a superposition of two Ca+ ions [TT]. In the 
experiment, the energy difference between the M = 1/2 
and 5/2 substates of the 3d ^£> 5/2 multiplet, monitored 
over 23 h served as a probe of Lorentz-violating effects. 

The anisotropic BBR shift produces a differential shift 
between M = 1/2 and 5/2 states mimicking the Lorentz- 
violating effects. Thus, anisotropic BBR is an important 
systematic effect. It was demonstrated in m that a 
factor of 10® higher sensitivity to Lorentz violation may 
be achieved with a similar experimental scheme with Yb"*" 
by monitoring the (4/^® 6s^) ‘^Ft/ 2 ^m= 7/2 - ^£ 7 / 2 .m=i /2 
frequency difference. Since this experiment will probe 
LLI at much higher sensitivity, study of anisotropic BBR 
is needed as it can be a major systematic effect for such 
an experiment. 

We note that singly ionized ytterbium ^’’^Yb+ with 
ultranarrow optical '^Si /2 - ^d? 3/2 and ^S' 1/2 - ^£ 7/2 transi¬ 
tions is also being pursued for a realization of an optical 
atomic clock and search for the temporal variation of the 
fine-structure constant a and the proton-to-electron mass 
ratio mp/me HIS]. 

The problem of anisotropic BBR shift of energy lev¬ 
els is practically unexplored so far, but may cause sys¬ 
tematic effects in a variety of experiments. In this 
work, we derived a general formula for the BBR shift 
of an energy level produced by a point-like source. A 
generalization to a finite source is obtained by integra¬ 
tion over angles of emitted thermal photons. The re¬ 
sult is expressed in terms of the scalar and tensor po¬ 


larizabilities of the atomic level. We also performed 
numeric calculation of the BBR frequency shifts for 
the Ca’*' ^£* 5 / 2 ,m= 5/2 " ^£* 5 / 2 ,m=i /2 transition and Yb'*' 
^£ 7 / 2 ,m= 7/2 “ ^£ 7 / 2 ,m=i /2 transitions due to their rele¬ 
vance to searches for Lorentz violation. 

An interaction of an atom in the state |0) with the 
electric field of a thermal photon emitted to a solid angle 
dfl leads to a blackbody radiation shift of an energy level. 
After integration over photon frequency, the BBR shift 
of the energy level |0) can be written as 


dE 

dVL 


3 

e i,k—l 


( 1 ) 


Here, we use a three-dimensional transverse gauge for 
photon polarization = (0, e) with polarization e nor¬ 
malized to the unit. Since photons are transverse, in this 
gauge ke = 0. The elements of the symmetric tensor aik 
are defined as 

_ (0Mi|TO)(m|4|0) , . 


where d = —r is the electric dipole moment operator and 
ujmo = Em — Eq is the difference between energy levels 
of the intermediate and |0) states. We use atomic units, 
i.e., |e| = = me = 1. An explicit form of the factor A 

is not important for the following derivation and we will 
restore it later. 

In the following we discuss only BBR effect caused by 
the electric field. The BBR caused by a magnetic field 
was considered in Ref. [TH] for a number of monovalent 
ions and proved to be negligible. Using the multipolar 
theory of blackbody radiation, developed in m , one can 
show that for the transitions in the Ca’*' and Yb+ ions, 
which will be discussed below, this effect can also be ne¬ 
glected. 

The electric dipole static polarizability of an atom in 
the state |0) is defined as 


Q^pol — ^zz 


2 E 


l(0|4|m)p 

^mO 


( 3 ) 


It can be conveniently decomposed into scalar and tensor 
parts: Opoi = as + at with the scalar polarizability 
given by 


as 



2y^ |(0|d|m)p 
3 ^ WmO 


( 4 ) 


The summation over photon polarizations in Eq. 0 is 
carried out using 
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where n = k/k. Then, Eq. Q is reduced to 


dE 

dQ 


i i,k 


= A 


3as ^ 


OCikTiiTl}- 


i,k 


The BBR shift in our case depends on the angle 9 be¬ 
tween the direction of the photon momentum k and the 
quantization axis z, defined by the direction of the mag¬ 
netic field. It is convenient to choose the vector k in 
the xz plane, i.e., ky = 0. Taking into account Eq. Q 
and noting that the product of the matrix elements 
(0|da;|n)(n|d^|0) = 0 and, respectively, a^x = ol^z = 0, 
we obtain 

A [30s O^xx^x ^zz'^z^ 

A [3q!s - ax:x sin^ 9 - a^z cos^ 9] . 


dE 


Taking into account that axx = cxyy and using the nor¬ 
malization condition + efy + = 1 and Eq. (j^ , after 

simple transformations, we obtain 



where 9i is the angle between the photon polarization 
vector £i and the z axis. 

Summing up over index i in Eq. (101, and using con¬ 
dition 


cos^0i -I- COS^02 + cos^0 = 1, (11) 

which is valid because the vectors ei, e 2 , and k are mu¬ 
tually orthogonal, we arrive at Eq. 


III. ANISOTROPIC BBR SHIFT FOR 

transitions 


Accounting for the fact that axx = ctyy and, hence. 


^XX 


3 q!s - azz 
2 


( 5 ) 


we express dE/dO, through Og, at, and cos^0. After sim¬ 
ple transformations we arrive at 


dE 

dVL 


= A' 



1 — 3cos^0 


4 


at , 


( 6 ) 


where A' = 2A. The factor A! can be easily determined, 
if we note that after integrating over d£l = sin0 d9 d(j) 
the second term in Eq. § disappears and we obtain 
AE = Ai:A!as. On the other hand, we have to arrive 
at the standard formula for isotropic BBR shift which, 
neglecting dynamic corrections, is given by m 

AE = -■^{aTT)^T^as, 


where the temperature T is given in a.u.. Finally, we 
obtain 


dE = - — 

15 

i-dcos^e 3M2-j(j-i-i) 

4 J(2J-1) 

Here we represent the tensor part at by 


-a2 


at — 


3M2 - J{J + \) 
J(2J- 1) 


-a2, 




( 8 ) 


where a^ is the tensor polarizability of the state |0). 

It may be instructive to present a different derivation 
of Eq. (|^, starting again from Eq. ([^. Assuming the 
polarization vectors ei ,2 to be real we can write 


dE 

dVL 


2 

(oLxxf'ix A ayyCty + ’ 

i=l 


( 9 ) 


We now apply Eq. ^ to the case of a \ JM) — \JM') 
transition between the ionic or atomic Zeeman sublevels. 
As we discussed above, the isotropic BBR shift, propor¬ 
tional to the scalar part of the polarizability, is very small 
for such a transition because it results only from a small 
difference between \JM) and \JM') energy levels. The 
main effect comes from the tensor part of the polarizabil¬ 
ity. 

Using Eqs. 0 and (|^ we write the \JM) — \JM') 
transition frequency BBR shift dEt as 


dEt = dEjM — dEjM- 
{anfT^ 
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(3cos"^0 — 1) 


m2 - M'^ 


dVL 


J(2J-1) 


Integration of Eq. (12) over fixed solid angle Hi leads 


to the BBR shift, corresponding to a maximal anisotropy 
100%, when the BBR is emitted to this solid angle and 
there is no BBR from the remainder. 

Let us consider a more realistic case when a certain 
portion of photons is emitted to the solid angle Hi at 
the temperature Ti and another portion of photons is 
emitted to the solid angle H 2 at the temperature T 2 , so 
that Hi -I- H 2 = 47r. The corresponding differential BBR 
shifts, which we designate as dE^^'^ (k = 1,2), are given 
by Eq. (pi) with T = 


Then, the total BBR shift can be found as 


AEt = 




1- 

H 

/o 

dH 

JQi 

dH 


'duA) 

dE[^^' 

dH. 


dH 

dH 


dfl 


(13) 


Performing integration in Eq. (13) over azimuthal an¬ 


gle (p from zero to 27r and over 9 from zero to a fixed 
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value 01 , we obtain 

(a7r)3(T4 - T^) 

' 20 ' J(2J- 1) 

X COS0i(l — cos^0i). 


AEt 


(14) 


As seen from Eq. (14), Ai?t turns to zero if Ti = T 2 , as 


it should be because it corresponds to the isotropic BBR 
case. 


It follows from Eq. (14) that the BBR shift is ^ 9f 


when 01 is small. Thus, this shift is greatly reduced with 
decrease in the solid angle to which the thermal photons 
are emitted. On the other hand, the anisotropic BBR 
shift is equal to zero when 0i = 90°. When 0i changes 
from 30 to 60°, cosdi(l — cos^0i) changes from 0.22 to 
0.38. 

Below we consider the anisotropic BBR shifts for the 
^A* 5 / 2 ,m= 5/2 - ^A* 5 / 2 ,m=i /2 transition in Ca’*' and for the 
^^ 7 / 2 ,M= 7/2 - ^^" 7 / 2 ,M=i /2 transition in Yb+. 


A. Anisotropic BBR shift for the 

^ 155 / 2 ,M=5/2 - ^l?5/2,M=i/2 transition in Ca+. 



FIG. 1: (Color onli ne) Dependence of AEt from the angle 
di for Ca”*" (see Eq. (15l) is represented by (blue) dash-dot- 
dotted line for Ti = 500 K, by dashed line for Ti = 420 K, 
and by dot-dashed line for T) = 350 K. Dependence of Ai?t 
from 61 for Yb+ (see Eq. (20l) is represented by (red) solid 
line for Ti = 450 K. The temperature T 2 = 300 K in all cases. 


In a recent paper m, the 


5/2,M=5/2 




b/2,M=l/2 


transition in Ca"*" was used to search for Lorentz invari¬ 
ance violation at a level comparable to the ratio between 
the electroweak and Planck energy scales. 

Using Eq. (14), we estimate the anisotropic BBR shift 
for this transition. The most accurate value of the ten¬ 
sor polarizability for the state was obtained in 

Ref. [T5], a 2 = —24.51(29) a.u.. To illustrate dependence 
of the BBR shift from the temperatures Ti and T 2 , we 
find AEt for three values of Ti (500, 420, and 350 K) 
and T 2 = 300 K. Substituting these values in Eq. ( [I4| , 
taking into account that for the room temperature 300 K 
(a 7 r)^T ‘*/20 « 4.9069 x 10“^® a.u., and expressing final 
results in Hz, we obtain 

AEt ~ cos 0 i(l — cos^ 0 i) 

-0.319 (Hz), Ti = 500K, 

-0.135 (Hz), Ti = 420 K, (15) 
-0.040 (Hz), Ti = 350 K. 


These dependences of AEt from the angle 0i in Ca®' 
^-D 5 / 2 m= 5/2 " ^^ 5 / 2 ,M=i /2 transition are illustrated in 
Eig.|^by three (blue) lines. As expected, AEt is equal to 
zero at 01 = 0 and 180°. The angle 0i = 180° corresponds 
to isotropic radiation. AEt also crosses zero when 0 i is 
90°. This case corresponds to isotropic radiation in the 
upper hemisphere. The energy shifts due to the BBR 
effect are largest for 0i « 55°. 


B. Anisotropic BBR shift for the 

^F’7/2,m=7/2 - ^E7/2,m=i/ 2 transition in Yb+. 

Recent work m identified several factors affecting 
the precision of the local Lorentz invariance tests with 


trapped ions. The two most important factors are the 
lifetime of the excited atomic state used in the Lorentz in¬ 
variance probe and sensitivity of this state to the Lorentz 
invariance violation effect, i.e., the size of the matrix el¬ 
ement of the corresponding operator. Both features are 
supplied by the metastable 4/^®6s® ^^ 7/2 state of the Yb+ 
ion, and the ^E 7 / 2 ,m= 7/2 - I 2 ,M=\I 2 transition was pro¬ 

posed as the probe of Lorentz-violating effects. To esti¬ 
mate this BBR shift we need to evaluate the value of the 
tensor polarizability 0:2 (^^ 7 / 2 ). We carried out calcula¬ 
tions in the framework of 15-electron configuration inter¬ 
action (Cl) method, following the approach described in 

Ref. [10]. 

The main features of this approach are briefly de¬ 
scribed below. All electrons are divided into the core 
and valence electrons. In our case [ls®,...,5p®] are the 
core electrons while 15 outer electrons belong to the va¬ 
lence subspace. 

In the framework of the Cl method we solve the eigen¬ 
value problem 

iLci$ = ^ci$, (16) 

where the many-electron wave functions $ belong to the 
valence subspace and are presented as a linear combina¬ 
tion of Slater determinants, 

$ = ^C/|det 7 ). (17) 

det/ 

The Cl Hamiltonian can be written as 

i7ci=i^core+ 51 H (18) 

N core j core 

where Ncore is the number of core electrons, E^ore is the 
energy of the core which includes kinetic energy of the 
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core electrons, Coulomb energy of their interaction with 
the nucleus and potential energy of the core-core electro¬ 
static interaction. The core-valence interaction, kinetic 
energy of the valence electrons and their interaction with 
the nucleus are included in the one-electron operators 
hP. The last term in Eq. (18) accounts for the interac¬ 


tion between valence electrons. 

We start from solving the Dirac-Fock equations for the 
[ls^,...,4/^^6s^] configuration. Then the 6pj and 5dj or¬ 
bitals are constructed for the 4/^^6s6p and Af^^6p5d con¬ 
figurations, correspondingly. The basis set used in the 
Cl calculations included also virtual orbitals up to 8s, 
8p, 7d, 7/, and 5g. We form configuration space by al¬ 
lowing single and double excitations for the odd-parity 
states from the 4/^‘*6p, 4/^^6s^ and 4/^^5d6s configu¬ 
rations and for the even-parity states from the 4/^‘^6s, 
4/^^6s6p and 4/^^6p5(i configurations to the orbitals of 
the basis set listed above. 

Solving the relativistic multiparticle Schrodinger equa- 
gives the eigenvector of the 4/^^6s^ 


7/2 


tion, Eq. ( p^ 

state whicnwe use to determine the tensor polarizability 
of this state. 

Using formalism of the reduced matrix elements we can 
write the expression for the tensor polarizability of the 
state $0 with total angular momentum J as [21] 


““ = '' (6(2J . 




J+J„ 


5J(2J-1) \ 

-^3)(2J+^l)(J-fl)j 
J 1 

1 J 2 f K-K ^ 


(19) 


where Jn is the total angular momentum of the interme¬ 
diate state 

A direct summation over all intermediate states 


in Eq. (19) requires a knowledge of the complete set of 


eigenstates of the Hamiltonian (16). Practically, this is 


impossible when dimension of a Cl space exceeds few 
thousand determinants, as in our case. To find the elec¬ 
tric dipole tensor polarizability of the 4/^^6s^ ‘^F 7/2 state 


we use the method of solution of an inhomogeneous equa¬ 
tion, described in detail in |21j . The random-phase- 
approximation corrections are also included. 

The result of our computation of the tensor polariz¬ 
ability is 0 f 2 (^F'ji 2 ) ~ —2 a.u.. Using this value and Ti = 
450 K and T 2 = 300 K, we obtain from Eq. (14) 


AEt 


—0.015 cos 6 * 1(1 — cos^6*i) Hz. 


( 20 ) 


The dependence of AEt on the angle 9i is shown in 
Fig. 0 by a (red) solid line. While the general behav¬ 
ior of AEt for Yb"*" is the same at 9i = 0, 90, and 180° 
as in Ca"*", the AEt is much smaller in the Yb'*' transi¬ 
tion considered here than in the Ca^ one. Suppression of 
the anisotropic blackbody radiation in the ^E 7 / 2 ,m= 7/2 - 
^^ 7 / 2 .M=i /2 Yb"*" transition is due to compactness of the 
Yb“'" 4/ orbital, resulting in the value of tensor polariz¬ 
ability, which is an order of magnitude smaller than that 
for the ^E 5/2 state of Ca'*'. Therefore, the anisotropic 
BBR shift is strongly suppressed for transition between 
substates of the Yb+ ‘^F'j /2 multiplet, and for Ti = 450 
and T 2 = 300 K its maximal (absolute) value at 6*1 « 55° 
is equal to 5.8 mHz. 

To conclude, we derived the formula for the anisotropic 
BBR shift of an energy level and performed numerical 
calculations of this effect for Ca’*' and Yb+ transitions of 
interest for study of Lorentz violation. We demonstrated 
that this effect strongly depends on the magnitude of 
the tensor polarizability of the level. In high-precision 
experiments, the anisotropic BBR can be a major 
systematic effect that should be specifically addressed in 
determining experimental uncertainties. 
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